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From first principles calculations we unravel surprisingly strong interactions between arsenic ions 
in iron-pnictides, the strength of which is controlled by the Fe-spin state. Reducing the Fe-magnetic 
moment, weakens the Fe-As bonding, and in turn, increases As-As interactions, causing giant re- 
duction in the c-axis. For CaFe2As2 system, this reduction is as large as 1.4 A. Since the large 
c-reduction has been recently observed only under high-pressure [l|, our results suggest that the iron 
magnetic moment should be present in Fe-pnictides at all times at ambient pressure. Finally, the 
conventional electron-phonon coupling in the collapsed phase of CaFe2As2 gives a maximum Tc of 
0.6 K and can not explain the ~ 12 K superconductivity observed in some experiments. Implications 
of these findings on the mechanism of superconductivity in iron-pnictides are discussed. 

PACS numbers: 74.25.Jb,67.30.hj,75.30.Fv,75.25.tz, 74.25. Kc 
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The recent discovery of superconductivity at Tc's up 
to 55 K in iron-pnictide systems [2, [sj, Ij, |5| has sparked 
enormous interest in this class of materials. So far two 
types of materials have been discovered. The first one is 
the rare-earth pnictide oxide layered systems, REOFeAs 
which is denoted as "llll";^ [3,, J, S, i]. The second 
class is the so called "122" systems with the chemical 
formula MFcaAsa (M=Ca,Sr, etc)[i, 0, i, i, 0. The 
122 systems are simpler in terms of their crystal structure 
due to absence of REO-layers. The crystal structure of 
CaFe2As2 is shown in the inset to Fig. 1. 

The pressure-induced superconductivity in these sys- 
tem is particularly interesting ll|, Il2l [l3| because it pro- 
vides a new avenue for investigation of the mechanism 
of the high-Tc superconductivity. Very recently a high- 
pressure neutron scattering study has reported that the 
CaFe2As2 system undergoes a surprising transition to a 
"collapsed tetragonal phase" (cT-phase) under appHed 
pressure in which the c-parameter is reduced from 11.7 
A to 10.6 A[l|. Motivated by this interesting report, in 
this letter we present a detailed first-principles study of 
Fe-pnictides with many surprising results. 

We discovered that the optimized c-lattice parameter 
of CaFe2As2 varies by about 1.5 A depending on the 
magnetic configuration considered. Such a giant cou- 
pling of spin-state of an ion with its lattice is almost 
unheard of and deserves detailed investigation. Reducing 
Fe-moment by about 1/2 gives lattice parameters that are 
quite close to high- pressure data[jj. Hence, the coUapsed- 
phase does not necessarily mean the total loss of Fe mag- 
netism as suggested in ref.[l| but could be due to large 
reduction of the Fe-moment. We explain this strange be- 
havior by noting that in iron-pnictide systems there are 
significant As-As hybridization (both intra- and inter- 
plane arsenic ions) whose strength is controlled by the Fe- 
spin state. Reducing the Fe-spin state reduces the Fe-As 
bonding, which in turn, increases the As-As bonding and 
causes the observed huge reduction in the lattice parame- 
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FIG. 1: (color online) Total energy as the Cal22 tetrag- 
onal phase goes to collapsed-tetragonal phase without any 
energy barrier during non-spin polarized structural optimiza- 
tion. The insets show the initial (T-Phase) and final (cT- 
phase) Cal22 structures with relevant bond-distances (in A) 
and angles (in degrees). 



ters. This effect is maximized in the case of Cal22 system 
due to close proximity of two arsenic ions in adjacent Fe- 
planes (see Fig. 1). We argue that since there is no ob- 
servation of large c-axis reduction during the normal to 
superconducting phase transition in iron-pnictides, the 
iron magnetic moment should be present at all times 
at ambient pressure. Otherwise one should have seen 
a large reduction in the c-axis as recently reported from 
high-pressure neutron scattering experiments [Ij. Hence 
it seems that the Fe-magnetic moment in its paramag- 
netic state (i.e. no long range order but non-zero Fe-site 
moment) is either required for the superconductivity or 



TABLE I: Various optimized structural parameters for NM, 
AF2, and API spin configurations, respectively. The exper- 
imental data are taken from Refs;&, il4i, 1^, [lg|. *Tlie AFl 
configuration goes to NM for Cao.5Nao.5Fe2As2 . 
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CaFe2As2 | 


NM 


5.63 


5.63 


10.39 


0.36251 


2.309 





0.0 


AFl 


5.65 


5.65 


10.60 


0.36440 


2.338 


1.3 


-16 


AF2 


5.61 


5.48 


11.61 


0.36695 


2.367 


2.2 


-100 


Exp. 


5.68 


5.68 


11.73 


0.3665 


2.370 


1.0 
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Cao.5Nao.5Fe2As2 | 


NM 


5.59 


5.59 


10.52 


0.36284 


2.31 





0.0 


AFl* 


5.59 


5.59 


10.52 


0.36284 


2.31 





0.0 


AF2 


5.43 


5.53 


12.05 


0.36536 


2.382 


2.4 


-97 


Exp. 


5.42 


5.42 


11.86 


- 


- 


0.0 
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BaFe2As2 | 


NM 


5.58 


5.58 


12.45 


0.3479 


2.319 





0.0 


AFl 


5.64 


5.64 


12.73 


0.35231 


2.382 


2.1 


-80 


AF2 


5.70 


5.59 


12.83 


0.3549 


2.408 


2.4 


-169 


Exp. 


5.52 


5.52 


13.02 


0.3545 


2.397 


1.0 
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LaOFeAs | 


NM 


5.64 


5.64 


8.59 


0.35944 


2.332 





0.0 


AFl 


5 .69 


5.69 


8.71 


0.35128 


2.393 


2.1 


-86 


AF2 


5.67 


5.73 


8.72 


0.34860 


2.409 


2.4 


-190 


Exp. 


5.70 


5.70 


8.737 


0.3479 


2.407 


0.35 


- 



at least is not detrimental to the superconductivity. 

The first-principles calculations were performed within 
the plane-wave implementation of the Perdew-Burke- 
Ernzerhof (PBE) generalized gradient approximation 
(GGA) to density functional theory as implemented in 
the PWSCF package [17[. The details of the parameters 
used in the calculations can be found in RefsllSl. lltA In 
order to get a general understanding of the iron-pnictides, 
in this study we consider one example of each class of 
pnictides; namely CaFe2As2 for the 122 system with 
the smallest Ca-ion available and the BaFe2As2 with 
the largest metal Ba. For the 1111 system, we con- 
sider LaOFeAs. We also study a doped 122-system, i.e. 
Nao.5Cao.5Fe2As2. Since in our ^/2 x V2-cell we have 
four chemical formula, we consider a supercell where two 
Na and two Ca are ordered. For each given system, we 
have performed full structural optimization including the 
lattice parameters and the atomic positions. We consider 
our optimization is converged when the maximum force 
on each atom is less than 0.005 eV/A and the pressure 
is less than 0.1 kbar. We have performed the full struc- 
tural optimization for non- magnetic (NM), i.e. "non-spin 
polarized", checkerboard antiferromagnetic (AFl) and 
stripe (AF2) spin configurations. Our results are sum- 
marized in Table 1 . As expected, the ground state for all 
four systems is the stripe AF2 phase and the optinrized 
parameters are in good agreement with the experimental 
data at ambient conditions. 

The most striking and surprising finding listed in Ta- 
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FIG. 2: (color online) The projected electronic density of 
states (PDOS) for As p-orbital (top-left) and Fe-d orbital 
(bottom left) for T- (black) and cT-phases (red). The right 
two panels show the band structure along F-M-F. The As- As 
Pz-bonding orbital is indicated by the arrow on the top left 
panel and it is shown in Fig. 3 



ble 1 is the giant dependence of the optimized c-lattice 
parameter on the spin-configuration considered. For the 
case of CaFe2As2, we note that AFl configuration is the 
next stable state (after the AF2) but the c-value is sig- 
nificantly reduced; 11.63 A versus 10.60 A for AF2 and 
AFl spin configurations, respectively. The difference is 
even larger, when the Fe-magnetism is ignored (i.e. non- 
spin polarized calculations). The optimized c-value for 
NM-state is 10.39 A, which is 1.34 A shorter than the 
experimental value at ambient pressure. We note that 
the optimized lattice parameters, a=5.65 A and 0=10.6 
A for the AFl phase are in reasonable agreement with 
the neutron data (a=5.8 A and c=10.6 A) 1]. Hence, 
based on this result, it is premature to conclude that the 
observed experimental cT-phase is due to the complete 
disappearance of the Fe-magnetism but as we shall see 
below it is closely related to the magnitude of the Fe- 
spin. 

In Fig. 1 we show that the T-phase directly collapses 
into the cT-phase without any energy barrier if the Fe- 
magnetic moment is ignored. During the c-axis collapse, 
there is significant and comparable decrease in the height 
of the Fe-As and As-Ca-As planes, indicating that the 
whole lattice almost uniformly shrinks. Since there is no 
energy barrier between the T-phase and the cT-phase, 
as soon as we loose the Fe-magnetic moment for some 
reason, we should see the expected c-lattice reduction. 
Since this does not seem to happen in the superconduct- 
ing samples, it is tempting to conclude that we have the 
Fe-moment present in the superconducting phase. This 
is a quite important result by itself to understand the 
mechanism of the high T^ superconductivity observed in 
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FIG. 3: (color online) Contour plot of the As- As p^- 
hybridization orbital for cT-phase (left) and T-phase (right), 
respectively. Note that the As- As hybridization present in 
both phases is much more significant in the cT-phase. There 
is also clear hybridization between intra- As atoms below and 
above the Fe-plane in the T-phase (right). 



iron pnictides. 

In order to get a better understanding of the cT-phase, 
we have performed electronic band structure calcula- 
tions and studied the projected atomic density of states. 
Our results for the T-phase agrees well with previous 
calculations 20|, |2l|, [22, l23[ . Fig. 2 shows that despite to 
the 1.4 A c-reduction, the band structure and the PDOS 
of both phases are surprisingly similar. The main differ- 
ence between the two phases is that most of the bands are 
shifted in lower energy in the cT-phase which is expected 
since the cT-phase has the lower energy. This downward 
shift of the bands is most obvious in the band structure 
plot along the T — M direction as shown in Fig. 2. In 
the T-phase, there were several states above the Fermi 
level which crosses the Yip along the M — T direction. 
This results in a larger density of states at the Fermi 
level. However for the case of cT-phase, we note that 
most of the bands just above E^? in the T-phase are now 
just below the E^? in the cT-phas and there is only one 
band which crosses the Yip along the V — M direction. 
This explains why N(£^f) is reduced in the cT-phase. 
The other and probably the most important difference 
is the change in the peak shape of the states near 3 eV 
below Ef. Visual inspection of these orbitals indicate 
that there is significant As- As hybridization in this sys- 
tem. The As- As hybridization becomes more significant 
in the cT-phase. This observation is quite unexpected 
and and as we shall see below it actually explains the 
mystery why the Fe-spin state has a huge effect on the 
As-z coordinates as first noted in 
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In order to demonstrate that there are large hybridiza- 
tion between As ions in the Cal22 system, we show the 
contour plots of the relevant orbitals in Fig. 3. It is clear 
that the As ion below the top Fe-plane makes a bond (or 
hybridizes) with the arsenic ion which is above the lower 
Fe-plane. Hence this overlap of the As- As along the c- 
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FIG. 4: (color online) Various structural parameters as a 
function of Fe-moment from fixed-spin calculations with fer- 
romagnetic iron spin configuration. 



axis makes this system quite isotropic and far from being 
layered system. According to bond-population analysis, 
the As-As bond strength increased almost twice in the 
cT-phase. Due to close proximity of the As ions in adja- 
cent Fe-layers, the observation of the As-As interaction 
is probably not that surprising. What is surprising is to 
see that there is almost the same type of hybridization 
between two arsenic ions on the same Fe-plane as shown 
in the right panel of Fig. 3. 

Since we have shown that the As ion above the Fe- 
plane has a strong overlap with the As ion below the same 
iron plane, their interaction is automatically increased as 
the Fe-As interaction decreases due to decrease in the 
Fe-moment which changes the chemistry of the Fe ion. 
Therefore, we have now a mechanism which explains why 
the As ion z-values get shorter with the decreasing Fe- 
moment. Our mechanism also explains why we see a 
smaller reduction in the c-axis for the LaOFeAs than the 
122 system as listed in Table 1. The reduction in the 
c-axis in the LaOFeAs system is due to the intra-plane 
As-As interaction only since there are no two adjacent 
FeAs planes to interact to with each other as in the case 
of Cal22. Our theory also predicts that for larger ions 
like Ba, we should see less c-reduction because the As-As 
distance between two adjacent planes are now larger due 
to larger ionic radius of Ba. In Table 1, we also show 
that similar c-reduction with Fe-spin occurs in the doped 
Nao.5Cao.5Fe2As2 system as well. 

In order to convince the reader further that the spin- 
state of the iron is the key parameter that controls the 
As-As bonding in these systems, we have performed fixed 
spin calculations for CaFe2As2 and the results are shown 
in Fig. 4. We consider ferro spin configuration to show 
that the spin-structure is a secondary effect and the main 
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FIG. 5: (color online) Phonon density of states (DOS), 
Eliashberg function {a'^F{uj)) and the frequency-dependent e- 
ph coupling A((jj) (dashed line) for CaFe2As2 in the cT-phase. 



effect is the on-site Fe-spin state. As we see from Fig. 4, 
when the Fe-moment is zero, the c-axis is the smallest 
consistent with the strongest As- As interaction (because 
As-Fe bonding is weak for non-magnetic Fe). As the 
Fe moment becomes significant, the Fe-As interaction 
gets stronger and therefore the As-As distance starts to 
increase to optimize the Fe-As bonding. Ironically at 
the Fe-moment of /xb=2.2, we have the ideal tetragonal 
arrangement of the four As ion around each iron dur- 
ing which the two As-Fe-As angles are equal. At this 
point, the low- to high-Fe spin transition is obvious dur- 
ing which the a-axis is reduced while c-axis is increased 
significantly. 

Since our results suggest that Fe-magnetism is either 
totally lost or reduced by half in the cT-phase, one won- 
ders if the observed '^ 12 K superconductivity in the 
vicinity of the collapse cT-phase of CaFe2As2ll3J can be 
explained by conventional electron-phonon (e-ph) cou- 
pling? In order to address this question, we have cal- 
culated phonon spectrum and Eliashberg function from 
linear response theory[17|. We used basically the same 
method and the equivalent parameters that are used in 
Refl25l for LaOFeAs. Our results are summarized in 
Fig. 5 and very similar to those for LaOFeAs. We ob- 
tained a value of electron-phonon coupling A — 0.23 
and the logarithmically average frequency ujiog — 218 K, 
which gives Tc = 0.6 K using the Allen-Dynes formula 
with fi* = (i.e. an upper bound for Tc). Hence, 
the mechanism of superconductivity in the cT-phase of 
CaFe2As2 is likely unconventional and probably it is the 
same as in the other iron-pnictides. This finding together 
with our results summarized in Table 1 suggest that in 
the cT-phase the Fe-magnetism is not totally lost with 
pressure and but it is partially reduced. Hence it seems 
that there is an optimum strength of Fe-spin state that 
is required for high-Tc superconductivity. 

In conclusion, we have revealed surprisingly strong As- 
As interactions in iron-pnictides. The strength of this 
interaction is controlled by the Fe-As chemical bond- 
ing. Reducing the Fe-moment, reduces the Fe-As bond- 
ing, which in turn increases the As-As interaction along 



the z-axis, causing arsenic atoms on opposite sides of 
Fe-square lattice to move towards each other. This ex- 
plains the high sensitivity of the z-atom positions and the 
large reduction of the c-axis with the loss of Fe-magnetic 
moment. We showed that under external pressure, the 
high Fe-spin AF2 structure (i.e. stripe phase) should 
transform to a new structure with low Fe-spin state and 
significantly reduced c-axis. We think that this could 
be the phase recently observed by high-pressure neutron 
scattering[l]. This is also consistent with our finding that 
the pressured induced ~ 12 K superconductivity can not 
be explained by the conventional e-ph coupling. The Fe- 
magnetism is still needed for an unconventional mecha- 
nism. Since at ambient pressure, we do not observe large 
c-axis drops in the superconducting samples, we conclude 
that the Fe-magnetic moment should be present at all 
times in these systems, at least in 122 materials such as 
CaFe2As2. The giant couphng of the on-site Fe-magnetic 
moment with the As-As bonding that we have discovered 
here may provide a mechanism for the superconductivity. 

The author acknowledges useful discussions with M. 
A. Green. 
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